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Abstract 

We discuss the interplay between collective flow and density profiles, de- 
scribing light cluster production in heavy ion collisions at very high energies. 
Calculations are performed within the coalescence model. We show how col- 
lective flow can explain some qualitative features of the measured deuteron 
spectra, provided a proper parametrization of the spatial dependence of the 
single particle phase space distribution is chosen. 

PACS number(s): 25.75 -q, 25.75 -Ld 

Keywords: relativistic heavy ion collisions, coalescence model, clusters, transverse 
flow, effective temperature. 



In the course of a relativistic heavy ion collision a hot and dense fireball is created 
in the interaction region. Due to the high internal pressure it expands and cools 
down, finally disintegrating into hadrons. The emergence of a collective flow can be 
considered as a signature that actually an extended piece of hot and dense matter is 
formed. Among the products of the reactions a few light nuclei and antinuclei have 
been observed(see, for instance, the Quark Matter '96 proceedings |L| and references 
therein), a very surprising fact for such high collision energies. A common scenario 
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employed to explain these observations is based on the coalescence model (for a review, 
see 0). A large amount of data on composite particle spectra is accumulated now at 
intermediate collision energies p3], where the multifragmentation of nuclei is a most 
striking phenomenon. Although some trends in the fragment spectra are similar to 
light cluster spectra in relativistic collisions, the mechanism of cluster production is 
in general very different in these two energy domains. In particular, the coalescence 
picture does not really work for intermediate mass fragments. 

In this letter we examine the effect of collective expansion on the final spectra 
of clusters, assuming that they are produced via coalescence. We also underline 
how their measurements can give an additional information on the latest stages of a 
relativistic nuclear collision. 

In the attempt of describing light cluster production in heavy ion reactions at very 
high energies, one encounters a somewhat subtle problem. Because these composite 
objects, typically d, d, t and 3 He, are very loosely bound, they can only be formed at 
the very late stage of the reaction. This is because the system is then quite dilute and 
interactions with the environment are therefore rare, preventing the formed clusters 
from breaking-up. On the other hand, it is known that light nuclei cannot be formed 
by scattering nucleons in free space, even when the process is gentle enough, simply 
as a consequence of energy-momentum conservation. The formation of a bound state 
requires the presence of a third body which carries away an amount of energy equal 
to the cluster binding energy. It is also clear that the system cannot be arbitrarily 
dilute and there must be a density around which the formation process is optimized. 

In the present study we assume that the production process is governed by two 
distinct factors. The early stages of the reaction and therefore the way particles are 
produced and emitted are parametrized via a many-body phase space distribution. 
This is what we call source funcion and it represents the probability that A nucleons 
are emitted at a given phase space point. The source function is taken at a sufficiently 
late time, such that the conditions previously discussed are fulfilled. The probability 
that these particles form a bound state, is taken as the overlap between the cluster 
and the /L-nucleon wave functions. This framework has been quite commonly used 
and it is well described in j3j. 
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Assuming that nucleons are emitted uncorrelated, one can factorize the generic 
A-body distribution as a product of single particle ones. Denoting a phase space 
point and the corresponding measure as 

(— # — # \ i ^ f ' % ^ Pi • -i a / -i \ 

r ii Pi) i dxi = ^27r) 3 ' % ~ ' ' ' ' ' 

we can write the phase space distribution of mass- A clusters as 

fA(r,p)= [[dx i f(x i )VA{xi,...,x A -,r,p). (2) 

^ i=i 

This formula expresses the fact that, among all A-particle states, represented by 
Y\f =1 f(xi), some can become bound with a probability Va- The integration goes 
over all phase space points, where particles are emitted. The formation probability 
Va is obtained by squaring a corresponding quantum-mechanical amplitude, as done 
in p|, [|. Below we adopt an approximation motivated by comparing the ranges of 
variation in phase space of two factors, the single particle distribution / and the 
formation probability Va- Obviously the first quantity has a much bigger range of 
variation than the second one, especially when considering a large and hot system. 
This allows us to set f, \ = r and pt = p/A, for % = 1, . . . , A, and writing the general 
formula for the phase space distribution of mass-A clusters in the form J7| 

f A (r,p) = [f(f,p/A)} A (3) 

This expression is the starting point of our subsequent analysis. 

Let us now specify the shape of the nucleon distribution in phase space. We assume 
that the system is in local thermal equilibrium, characterized by a temperature To, 
considered to be constant throughout the whole fireball at the freeze-out stage of 
the reaction. We also assume that particles are subject to a collective velocity field, 
often also named collective flow. At very high energies it is generated by the partial 
transparency of nuclei, along the longitudinal (beam) direction, and by the pressure 
created in the hot overlap zone, in the transverse direction. Since the dynamics in 
these two directions is very different, we disregard possible correlations and represent 
the collective velocity field as a sum of two independent contributions, 

v(f) = v L (f L ) + v T {r T ) , (4) 
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where = z e z and tt — x e x + y e y . The nucleoli momenta in a local frame k obey 
a thermal distribution with temperature T . The transformation to a global frame is 
made with a boost of velocity v(r). It is well known that the longitudinal dynamics is 
highly relativistic, while the transverse expansion, in the first approximation, can be 
considered non- relativistic, at least for nucleons. Therefore, the nucleon momentum 
in the global frame can be written as the sum of thermal and flow components as 
Pt = £>r + m-yV(rr). In the following discussion we ignore all issues related to the 
longitudinal dynamics, focusing attention on the transverse direction. The transverse 
velocity field is parametrized as 

v T (r T ) = v f (jQ e T (5) 

where Vf and Ro are the strength and scale parameters of flow and the power-law 
profile is characterized by the exponent a. In building the phase space distribution, 
we follow ref. |J. Assuming cylindrical symmetry we represent the nucleon density 
in a factorized form 

p{r) = N n L (f L )n T {r T ) , (6) 

such that Hi and ut are normalized to 1 in the respective domains. Using (|3|) we can 
now calculate the cluster phase space distribution function. Because of (H) and (Bp, 
it also factorizes into longitudinal and transverse parts, namely 

fA(r,p) = f L (r L J L ) f T (r T J T ) . (7) 

The transverse contribution for clusters of mass number A is therefore given by 

(p T -Mv T ) 2 



f T (r T ,p T ) = (2nf^- )e - ™ Ta nA[ f Th (8) 

where M = Am is the cluster mass and 

n A (f T ) = N T (A) [n T (f T )} A (9) 

is the transverse part of the cluster density, with the normalization factor Nt(A). 
Position and momentum of particles, completely uncorrelated in a purely thermal 
system, are now partially linked due to the presence of collective flow. 
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The transverse momentum spectrum of clusters is obtained by integrating expres- 
sion (0) over the whole volume and around a particular value p L of the longitudinal 
momentum^. The pr-spectrum for clusters of mass number A can be written in the 
form 

dN A 
dp\ 

where va(j>l) ls the total number of clusters of mass A produced at p L and 



p l = va(Pl) Sa(pt) (10) 



is the pr-dependent part of the momentum spectrum. This last factor is quite in- 
teresting. If flow were absent, the integral would give directly the Boltzmann factor, 
but the present case is, in general, more complicated, and a numerical treatment is 
needed. The most common parametrization used in the literature combines a gaus- 
sian profile for the nucleon density || with a linear profile (a = 1) for collective 
flow (see [f|, |l0j, especially in relation with source parametrizations in interferometry 
studies). Only this choice allows for an analytical solution, which is the Boltzmann 
distribution 
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Pt 



but now with the modified effective temperature (slope parameter) 

T* = T + mvj . (13) 

At first sight, this result looks appealing, but it actually contradicts both intuition 
and experiment. What is wrong in the previous expression is the dependence of the 
slope T* only on m but not on M, as one would expect also by looking at the slopes 
extracted from measured spectra fTI| . When performing the integral in (ITTJ), one 
notices an interesting feature. From @ one sees that the density of clusters of mass 
A is proportional to the A-th power of the nucleon density. In the case of a gaussian 



a relativistic formulation, a more familiar notation in terms of rapidity y instead of pl would 
appear, without anyway affecting our discussion on transverse spectra. 
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profile, one can see that the A-cluster density shrinks towards the central region. 
This is easy to understand, since it is clearly more probable to make a cluster where 
there are many particles than on the tail, where there are only a few. Together with 
this, we choose a linear flow profile. This is parametrized in (f|) defining Vf as the 
flow strength at the surface of the density distribution, characterized by the scale 
parameter Rq. What happens with the gaussian profile is that the actual size of the 
cluster density has a smaller radius, thereby picking up a smaller value for the flow 
velocity at the surface, as compared to the case of single nucleons. This effect exactly 
cancels the A-dependence of M in (0). The other extreme would be to take a uniform 
density with a sharp surface at a given radius. Any power of this function would give 
the same profile, with the same radius. In other words it is equally probable to have 
clusters everywhere in the region with non-zero density. As a consequence we expect 
in this case that the slope parameter will depend on M, since the flow velocity at the 
surface is the same for all clusters. This is not the whole story. The flow profile could 
have a smaller exponent (a = 1/2, for example). Indications of such a behaviour have 
been observed in microscopic models of heavy ion collisions such as RQMD JTJ], 
Now some dependence of T* on M would appear, even for a gaussian density profile. 

Let us now look more closely at the interplay between flow and density profiles. 
It is clear that they cannot be considered independently because the density shape 
at a given time during expansion is the result of the particle motion characterized 
by the collective velocity field. The information about the profiles of density and 
collective velocity, can be extracted, in principle, from the energy spectra of different 
clusters 0. Unfortunately this is not an easy task because of the sensitivity of energy 



spectra to all kinds of corrections [14, 15|. We prefer a more global analysis where 



the effective temperature (slope parameter) is extracted from the A-dependence of 
the mean transverse energy[| 

<^>WW^§^,«. (14) 
For a classical Boltzmann gas at temperature Tq we have < Et >a = To, where the 



2 This approach was further developed in |l6| , including in the analysis also the variance of the 
energy distribution. 
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usual factor 3/2 has changed to 2/2 since we consider only the transverse degrees of 
freedom. In the present case we have instead 



<E T > A = <E T >* + <E T > A l , (15) 

where the first term corresponds to the purely thermal, Boltzmann gas, while the 
second contribution arises due to the presence of flow (It vanishes if we set Vf = 0). 
We define the effective temperature via 

T„ = < E T > A = T + AL n] (A) Mvj . (16) 

The coefficient in front of the flow term depends on the flow parameter a, is a func- 
tional of the transverse density and is a function of the cluster mass number A. 
Different choices of density and flow profiles will result in a different A-dependence. 
The explicit expression for this coefficient is 

A L™ ] (^) = \J d 2 r T (l^MrT) • (17) 

It can be calculated analytically for the two interesting cases of gaussian and box 
profiles for the density, for all values of a. In the first case one obtains 

A^ SS (A) = ^ T(a + 1) , (18) 

where T is Euler's Gamma function. One can readily see that for a linear flow profile 
(a = 1) the coefficient is equal to 1/A and it exactly cancels the A factor carried 
by M in (0). For lower powers of a the situation changes one maintains a weak 
A-dependence. In the case of a box profile we obtain 

A*™ (A) = , (19) 

a v ; 2a + 2 ' v ; 

independent of A, as we expected after our previous discussion. Choosing parame- 
ters according to Table [L], we illustrate the results in Figure [I], where the effective 
temperature is plotted as a function of mass number A. The higher curves represent 
the extreme case of a box-shaped density which gives the strongest A-dependence of 
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T*. The other extreme, as pointed out previously, is the gaussian density with linear 
flow, which gives an A-independent T*. Figure [TJ suggests that the choice of a gaussian 
profile for the density and a flow profile with a = 1/2 give the best agreement with 
the measured values of the slopes. 

There is another interesting feature regarding cluster spectra, which can be experi- 
mentally measured. In the early days of heavy ion physics the proportionality relation 
between cluster spectra and the corresponding powers of single particle spectra was 
quite well established 0. In recent experiments at much higher energies a momen- 
tum dependence in the proportionality constant Ba was observed [|17j]. Namely, Ba 
increases with increasing transverse momentum. In the present analysis we compare 
the transverse momentum spectrum of clusters of mass number A with the A th power 
of the single particle spectrum, 

dN A I ._ . / dN I \ A 

^h L = bA ^ PT) [d(^h)^ (20) 

where the first factor does not coincide with the usually quoted Ba because it is 
calculated for a small window around p L . We therefore discuss the pr-dependence in 
this factor. Again, we perform the calculations with box and gaussian profiles, taking 
a = 1, 1/2. Using fllCf) we obtain 

b A (p L ,PT)=C A (p L ) A , (21) 

[Sa{pt/A)\ 

where ca is a normalization factor which gives the order of magnitude of 6^, but 
does not affect its pr-dependence. In Figure |]we show various plots of b 2 , using the 
parameters from Table [3]. Although microscopic simulations are able to more or less 
reproduce this feature Ul8| , it is instructive to understand how it arises within the 
simple picture presented above. This behaviour is a pure manifestation of collective 
flow, which only cluster measurements can reveal. This effect depends on the relation 
between flow and density as we discussed above, resulting in turn in different shapes 
of momentum spectra for clusters and single nucleons. We emphasize again that 
the linear velocity profile and the gaussian shape for the density distribution are 
in contradiction with the pr-dependence of 6^ (In this case both cluster and single 
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particle spectra have the same slope). Also the choice a = 1/2 does not help much, 
suggesting that a better understanding of the density shape is necessary. Therefore 
we indicate in Figure ^| that surface formation of clusters, at a slightly earlier time 
with respect to the complete disintegration of the system, could improve ou scenario. 
This is done by performing the spatial integration over a spherical shell from Rq/2 to 
Ro and is equivalent to having a density with a depleted central region, as suggested 
in [ i~9fl . The actual situation is clearly a combination of this early surface emission 
and final bulk disintegration and a consistent implementation of this aspect, together 
with a proper description of time evolution, is the subject of our current study. 

In summary, we have shown that a suitable implementation of collective flow 
can account for important qualitative features of light cluster spectra, measured in 
heavy ion collisions at very high energies, even though more has to be done to build 
a consistent and quantitative description of the late expansion stage. The observed 
^-dependence of the slope parameters and the py-dependence of the coalescence coef- 
ficients impose serious constraints on the spatial profiles of the collective velocity and 
the particle density at the freeze-out stage. The most common parametrizations for 
both flow and density profiles fail to reproduce these features. Quantitative conclu- 
sions will be possible in the near future when cluster spectra for large and symmetric 
collision systems will be available. 
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Table 1: Parameters used to perform the calculations. They are not chosen in order 
to fit the data but only to show the qualitative beaviour in the following figures. 



Density 


a 


v f 


To(MeV) 


#o(fm) 


Box 


1/2 


0.63 


140 


8 




1 


0.72 


140 


8 


Gauss 


1/2 


0.48 


120 


8 




1 


0.34 


140 


8 
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Figure 1: Effective temperature T* as a function of mass number A. The top curve 
for each choice of density profile corresponds to a — 1/2, while the lower is for a — 1. 
The straight lines could perfectly coincide if a more accurate choice of Vf were made. 
This is evident from (16), due to the A-independence of for the box profile. 
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Figure 2: ft 2 factor as a function of transverse momentum. The top curve for each 
choice of density profile corresponds to a = 1/2, while the lower is for a = 1. The 
curve labelled "Surface" corresponds to integration over a spherical shell from Rq/2 
to Rq, in order to simulate surface emission. 
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